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The prospects of forming ultracold molecules in 2 S states 
by magnetoassociation of alkali-metal atoms with Yb 
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Durham University, South Road, Durham, DH1 3LE, United Kingdom 
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We explore the feasibility of producing ultracold diatomic molecules with nonzero electric and 
magnetic dipole moments by magnetically associating two atoms, one with zero electron spin and 
one with nonzero spin. Feshbach resonances arise through the dependence of the hyperfine coupling 
on internuclear distance. We survey the Feshbach resonances in diatomic systems combining the nine 
stable alkali-metal isotopes with those of Yb, focussing on the illustrative examples of RbYb and 
CsYb. We show that the resonance widths may expressed as a product of physically comprehensible 
terms in the framework of Fermi's Golden Rule. The resonance widths depend strongly on the 
background scattering length, which may be adjusted by selecting the Yb isotope, and on the 
hyperfine coupling constant and the magnetic field. In favorable cases the resonances may be over 
100 mG wide. 



I. INTRODUCTION 

The successes of cooling gases of atoms to ultracold 
temperatures have led to great interest in producing 
molecules at similar temperatures. Because molecules 
have a richer internal structure and more complex inter- 
actions than atoms, ultracold (/xK) molecules offer the 
possibility of exploring a wide range of new research ar- 
eas, including high-precision measurement [IH3], quan- 
tum information [5] [5] and quantum simulation [5] . 

Molecules may be formed in ultracold atomic gases ei- 
ther by photoassociation [7] or by magnetoassociation [8] . 
In the latter, cold atomic clouds are subjected to time- 
dependent magnetic fields that convert atom pairs into 
molecules by adiabatic passage across zero-energy Fesh- 
bach resonances [9]. Recent years have seen substan- 
tial progress in producing ultracold molecules made up 
of pairs of alkali-metal atoms [TUHPE] . The molecules are 
left in high vibrational states and are susceptible to col- 
lisional trap loss. For KRb [15], Cs2 [E], an d triplet Rb2 
[T7] , it has been possible to transfer the molecules to the 
absolute ground state by Stimulated Raman Adiabatic 
Passage (STIRAP). 

There is now great interest in the formation of cold 
molecules that have both electric and magnetic dipole 
moments [19-24 . Such molecules offer additional possi- 
bilities for manipulation, trapping, and control because 
they can be influenced by both electric and magnetic 
fields. In the present paper we investigate the prospects 
for magnetoassociation of alkali-metal atoms (Alk) with 
1 S atoms (specifically Yb) to form 2 E heteronuclear di- 
atoms with electron spin S = 1/2. 

Ytterbium is an excellent candidate for pairing with 
the alkali metals. It has 7 stable isotopes (5 zero-spin 
bosons, 2 fermions), and a closed-shell, singlet-spin elec- 
tronic structure. Both bosonic 25-281 and fermionic 
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[291 150] isotopes have been cooled to quantum degen- 
eracy. Different isotopic combinations have different 
scattering lengths, and produce molecules with different 
binding energies; they thus have Feshbach resonances at 
different magnetic fields. 

The existence of magnetically tunable Feshbach reso- 
nances requires coupling between a continuum scattering 
state of the atomic pair and a molecular state that crosses 
it as a function of magnetic field. For pairs of alkali-metal 
atoms, this coupling is provided by the difference be- 
tween the singlet and triplet potential curves and by the 
magnetic dipolar interaction between the electron spins. 
However, neither of these effects exists in systems of the 
type considered here. Instead, the most significant coup- 
ling between the atomic and molecular states is provided 
by the independence in the hyperfine coupling constant 
of the alkali- metal atom [19] . Such independences exist 
in alkali dimers |31| , but in that case they merely produce 
small shifts in bound state energies and resonance posi- 
tions, rather than driving new resonances. If the closed- 
shell atom has non-zero nuclear spin, it can also couple 
to the unpaired electron spin. For the case of LiYb [21] , 
this coupling has been found to be much stronger than 
that due to the Li nucleus. However, this latter effect is 
less important for the heavier alkali-metal atoms consid- 
ered here, where the coupling to the alkali-metal nucleus 
itself is stronger. 

In previous work, we extracted resonance positions and 
widths for RbSr [19] and LiYb [24] from coupled-channel 
quantum scattering calculations. In the present paper, 
we extend these studies to a range of heavier systems and 
show how the widths may be broken down into their con- 
tributing factors within the framework of Fermi's Golden 
Rule. 

The theoretical development presented here is applica- 
ble to any system made up of an alkali-metal atom paired 
with a closed-shell atom. In the present study, we have 
considered the whole range of Alk-Yb systems, but we 
focus our presentation on the illustrative examples of Rb- 
Yb, for which the scattering lengths are approximately 



known, and Cs-Yb, for which they are as yet unknown. 
In section Ql] we describe the theoretical methods used. 
In section |HI| we present our results, with discussion of 
system characteristics that lead to Feshbach resonances 
suitable for molecule formation. 



II. THEORY 

A. Collisions between alkali- metal and closed-shell 

atoms 

The Hamiltonian for an alkali-metal atom a in a 2 S 
state, interacting with a closed-shell atom b in a *S state, 

is 



H 



2// 



d 2 L 2 
AR 2 + R 2 



U(R) + H a + H b (1) 



where L is the two-atom rotational angular momentum 
operator and U{R) is the interaction operator. H a and 
Hb are the single-atom hamiltonians, 



H a = Caia ' S + I <3W'\/„. 

Hb = gb/J-Nib,zB, 



ge^BSz ) B (2) 

(3) 



where s, i a and ib are the electron and nuclear spin opera- 
tors, g e , g a and gb are the electronic and nuclear g factors, 
and /ib and /ijq are the Bohr and nuclear magnetons. Ca is 
the hyperfine coupling constant for the alkali-metal atom 
and B is the external magnetic field, whose direction de- 
fines the z axis. In the present work we use lower-case 
angular momentum operators and quantum numbers for 
individual atoms and upper-case for the corresponding 
molecular quantities. 

The interaction of a 2 S atom with a : S atom produces 
only one molecular electronic state, of 2 S symmetry. 
However, the hyperfine coupling constant of the alkali- 
metal atom is modified by the presence of the closed- 
shell atom [TH] , and if % ^ then there is also hyperfine 
coupling involving the nucleus of atom b [23] , 



Ca(R) = Ca + ACa(R); 

UR) = A( b (R). 
The interaction operator U(R) is thus 



(4) 
(5) 



U(R) = V(R) + A( a (R)i a ■ s + A( b (R)ib • S, (6) 

where V(R) is the electronic interaction potential. Most 
of the theory presented here remains applicable when 
atom a is a non-alkali-metal atom in a multiplet-S state. 
Figure [I] shows the energy levels of the 133 Cs atom, 
with i a = 7/2 (black solid lines). At zero field the lev- 
els may be labeled by quantum numbers f a ,mf >a , where 
f a = i a ± 1/2, whereas at high field the nearly good 
quantum numbers are m sa and m ia . In the present 
paper we indicate the lower and upper states for each 
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FIG. 1. (color online). Hyperfine energy levels for 133 Cs in 
its ground 2 S state (solid black lines) and for near-threshold 
states of CsYb arising from the upper hyperfine manifold 
|o!2, rrif,cs) (dashed blue lines), calculated for a potential with 
a scattering length of —38 bohr. The n — — 1 level is almost 
hidden by the threshold. The solid circles mark the Feshbach 
resonances that occur at crossings between bound states and 
atomic thresholds with the same value of m/,c s - 



mf, a = 

tively. 



>a + m it a as |ai,TO/,a) and |a 2 ,TO/, a ) respec- 



The Hamiltonian (TTj) may be written as the sum of a 

zeroth-order term H° and a perturbation H' , 



ff° 



tr 



d 2 

di? 2 



2// 
H' = A( a (R)i a 



L 2 
R 2 



V(R) +H a + H b - (7) 



s + A( b {R)i b ■ s. 



(8) 



The zeroth-order Hamiltonian is separable, and 
its eigenfunctions are products of atomic functions 
\on,mf ta )\ib,mi,b) and radial functions ip(R). The latter 
are eigenfunctions of the 1-dimensional Hamiltonian 



a 2 

2// 



d 2 
dR 2 



L(L + 1) 
R? 



V(R), 



(9) 



with eigenvalues E n . The eigenvalues of H° are E n + 
E a + Eb, where E a and Eb are the eigenvalues of H a 
and Hb- By contrast with the alkali-metal dimers, the 
molecular states thus lie almost parallel to the atomic 
states as a function of magnetic field. They also have al- 
most exactly the same spin character. The only terms 
in the Hamiltonian (fTl) that couple \ax,mf >a )\ib,mi t b) 
and \a2,m'f a \ib,m' ib ))sxe the weak couplings involving 
ACa(R) and A£&(i?). The former couples states with 



l f,a = m f,a 



and 



•i,b 



TTH,b: 



while the latter couples 



states with m't a = m/ >a ± 1 and m' ib = m i,b T 1- 

The Hamiltonian (IT]) is entirely diagonal in L, so reso- 
nances in s-wave scattering can be caused only by L = 



bound states. The only interactions off-diagonal in L are 
spin-rotation and nuclear quadrupole interactions, which 
are neglected in the present work. This again contrasts 
with the alkali-metal dimers, where the magnetic dipolar 
interaction between the electron spins and second-order 
spin-orbit coupling provide relatively strong interactions 
that produce resonances from bound states with L > 
in s-wave scattering. 

Figure [T] shows the highest few vibrational states for 
CsYb with spin character \a2,Tn'r a ) and vibrational 
quantum numbers n — — 1, —2 ... — 7 (with respect to 
threshold) as dashed blue lines, calculated for a poten- 
tial V(R) with an s-wave scattering length a = —38 bohr. 
The couplings involving A£ a (R) give rise to Feshbach res- 
onances at fields where bound states |«2, inf >a ,n) cross 
thresholds |ai,m/, a ), shown as solid circles in Figure [I] 
In the present work we neglect couplings due to A(b(R) 
and set % = for all isotopes. This will give accurate 
results for resonances with m'r a = mf >a but will suppress 



resonances with 
171 Yb and 173 Yb. 



l f,a 



m/ ja ±l, which actually exist for 
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FIG. 2. (color online). Electronic potential energy curves 
V{R) from CCSD(T) calculations on the Alk-Yb systems. 



System R e V(R e ) C 6 

(bohr) (mE h ) (E h qg) 

LiYb 6.65 -7.48 1594 

NaYb 7.61 -4.61 1690 

KYb 8.88 -3.36 2580 

RbYb 9.28 -2.99 2830 

CsYb 9.72 -2.83 3370 



B. Electronic Structure Calculations 

1. Potential Energy Curves 

We have constructed the electronic potential energy 
curves V(R) by carrying out electronic structure calcula- 
tions at short and medium range and switching to a form 
incorporating dispersion interactions at long range. 

We obtained ground-state potential curves for NaYb, 
KYb, RbYb and CsYb from CCSD(T) calculations 
(coupled-cluster with single, double, and non-iterative 
triple excitations) using the Molpro package [32] • For 
Yb, we used the quasi-relativistic effective core poten- 
tial (ECP) of Dolg et al. and its corresponding basis set 
[33j , with 60 electrons in the inner 4 shells represented by 
the ECP and the remaining 10 electrons (p 6 s 4 ) treated 
explicitly. ECPs pH] and their corresponding basis sets 
[55] were also used for K, Rb, and Cs. An ECP was not 
used for Na; all 11 electrons were represented with the 
cc-pvqz basis set of Prascher et al. 36 . For each sys- 
tem, CCSD(T) calculations were carried out at a series 
of points from 2 to 40 bohr and the potentials were then 
interpolated using the reproducing kernel Hilbert space 
(RKHS) method [37]. The resulting potential curves are 
shown in Figure[2] together with the LiYb curve of Zhang 
et al. |38j , obtained using similar methods but with a fully 
relativistic ECP for Yb. The well depths and equilibrium 
distances are given in Table [I] 

At long range the potential curves were represented as 



V \R) — — C&R — C&R — CiqR 



(10) 



TABLE I. Properties of the interaction potentials used in the 
present work. The well depths and equilibrium distances are 
from CCSD(T) calculations and the Cq coefficients are from 
Eq. Jill. 



Slater-Kirkwood formula, 



a 



ab _ Cl a Cl b a a (0)a\0) 

Q a (a 6 (0)) 2 + Cf(a°(0)) 5 



(11) 



The C@ coefficients used for the long-range potential were 
obtained from Tang's combination rule |39j based on the 



using the homonuclear Cq coefficients for Alk-Alk [4"D] 
and Yb-Yb [41] and the static polarizabilities a(0) for the 
alkali- metal atoms [ID] and Yb [52]. Equation (11) gives 
Ce coefficients well within 1% of the values of ref. [ID] 
for all the mixed alkali-metal pairs. The results for the 
Alk-Yb systems are included in Table [I] The Cs and Cio 
terms were omitted except when fitting to the experimen- 
tal spectra for RbYb as described in section |HI A| below. 
The short-range and long-range regions of the potential 
were joined using the switching function of Janssen et al. 
1551 between the distances 28 and 38 bohr. 



2. Hyperfine Coupling 

The hyperfine coupling constant of an atom is a mea- 
sure of the interaction between its nuclear spin and the 
electron spin density at the nucleus, which in the case of 
an alkali metal comes principally from the single valence 
electron. Approach of another atom perturbs the elec- 
tronic wavefunction and alters the spin density at the 
nucleus, so that the coupling between the electron and 
nuclear spins becomes a function of internuclear distance 
R. 
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FIG. 3. (color online). Distance dependence A£ a (-R) of the 
hyperfine coupling constants for the Alk-Yb systems. 



Co 


(MHz) /3 (bohr" 2 ) R c 


(bohr) 


6 Li 


-48.8 


0.0535 


4.92 


7 Li 


-129. 


0.053^] 


4.92 


23 Na 


-258. 


0.0553 


5.18 


39 K 


-42.4 


0.0474 


6.09 


40 K 


52.5 


0.0474 


6.09 


41 R 


-23.3 


0.0474 


6.09 


85 Rb 


-177. 


0.0357 


5.71 


87 Rb 


-597. 


0.0357 


5.71 


133 Cs 


-377. 


0.0260 


5.54 



a The value of /3 for Li was reported incorrectly in ref. 1241 . 

TABLE II. Parameters of the Gaussian functions used to 
represent A£ a (R), the distance-dependence of the hyperfine 
coupling constant of an alkali-metal atom interacting with 
Yb. 



We have calculated the hyperfine coupling constants 
Ca(R) for the Alk-Yb systems, using density-functional 
theory with the KT2 functional [44] . as implemented in 
the ADF suite of programs [45] . We fitted these results 
to a variety of functional forms and found that, in the 
range of R for which the vibrational wavefunctions are 
non-zero, a Gaussian function A( a (R) — Q e~^ R ~ Rc ^ 
gave an adequate fit to the DFT results. These functions 
are shown in Figure [3] for each of the Alk-Yb systems and 
the parameters are given in Tabic In] 



netoassociation is experimentally feasible. Defining B ZCIO 
as the field where a{B) = near resonance, Eq. ( 12 1 im- 
plies A = B zcro - B ICS . 

In the present work, we obtained scattering lengths 
a(B) principally from coupled-channel calculations. The 
coupled equations for each field B were constructed in 
an uncoupled basis set |s a TO Sja )|z Q r7ii ]a )|LMi) and solved 
using the MOLSCAT package g7][IB]. The s-wave scat- 
tering length was then obtained from the identity [H] 
a = (ifc) _1 (l — 5oo)/(l + <Soo)> where Sqq is the di- 
agonal S-matrix element in the incoming channel, k = 
h~ 1 (2/j,E co i) 1 ^ 2 , and the collision energy E co \ was taken 
to be 1 nK xfce- MOLSCAT has an option to converge 
numerically on the fields corresponding to both poles and 
zeroes in a(B), allowing the extraction of A. 

It should be noted that Eq. ( 12 1 characterizes the scat- 



tering length near resonance only for purely elastic scat- 
tering. If there exist lower-energy channels that allow de- 
cay, then a{B) has a non-zero imaginary component and 
does not follow the simple pole formula (12) [49]. For 
ib = 0, the Hamiltonian (TIT) allows only elastic scattering 
even when the alkali-metal atom is in a magnetically ex- 
cited state. However, when if, ^ 0, couplings involving 
ib ■ s can change rrif ta , and for alkali-metal atoms in mag- 
netically excited states this provides additional couplings 
to lower-lying thresholds. We have previously described 
the behavior of a(B) for resonances in such states for the 
LiYb systems [24] . 



D. Resonance widths from Golden Rule 

Coupled-channel calculations of resonance widths are 
straightforward but provide relatively little insight into 
the factors that affect resonance widths. We therefore 
develop here an alternative approach based on Fermi's 
Golden Rule that allows us to understand the factors 
that determine the widths. 

Fermi's Golden Rule gives an expression for the width 
of a Feshbach resonance in terms of the off-diagonal ma- 
trix element of H' (Eq. M) between the bound state 
\cK2,mf,a,n) (with vibrational quantum number n) and 
the continuum state \ai,mf^ ai k) (labeled by wavevector 
fc, where E co \ = h 2 k 2 /2/i). The Breit-Wigner width in 
the energy domain, Te, is 



Resonance widths from coupled-channel 
calculations 



Near resonance, the s-wave scattering length a{B) as 
a function of magnetic field B behaves as [115] 



a{B) = dbg I 1 



B - B rc 



(12) 



where B TOS is the resonance position and <Zb g is the back- 
ground scattering length. The magnitude of the reso- 
nance width, A, is critical for determining whether mag- 



T E (k) = 2ir 



&2,mf,a, 



H' 



Q-i,m f:a 



k 



(13) 



where the continuum function is normalized to a 
(5-function of energy and has asymptotic amplitude 
(2/j,/nh 2 k) 1 / 2 . At limitingly low collision energy, Tsik) 
behaves as [9], 



T E (k) -^-> 2fca bg r 



(14) 



where o,b g is the same background scattering length as in 
Eq. (12 1. Tq is independent of energy and is related to 



the magnetic resonance width A of Eq. ( 12 ) by 



A = 



5(j, T 



(15) 



where Sfi rcs is the difference between the magnetic mo- 
ment of the molecular bound state and that of the free 
atom pair, which is simply the difference in slope of the 
crossing lines in Figure IT] 

The expression for the magnetic resonance width A 
factorizes into spin-dependent and radial terms, 



where 



A = 



nl, 



fcabgO/ircs 



(16) 
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FIG. 4. (color online). Off-diagonal matrix elements of i a • $ 
for 133 Cs between pairs of hyperfine states with the same value 
of m/, . 



I mja (B) = (a2,mf,a\i a • s\a 1 ,m f . a ) (17) 

and 

i n k= MR)&UR)MR)dR. (18) 

Jo 

The quantity I m , a (B) is a purely atomic property, 
which arises because states |aj,m/ i0 ) are eigenfunctions 
of H a . Pairs of states with the same m/. a are coupled 
through the operator A( a (R)i a ■ s. At zero field, the 
states are eigenfunctions of i a ■ s, so that the perturba- 
tion has no off-diagonal matrix elements. At sufficiently 
high field, however, the states are well described by quan- 



tum numbers m s 

and s a — 1/2, 



and 



such that for a given m 



f,a 



-&- 



-&- 




in^w ^ — 



R (Bohr) 



Irn f JB) 



^ 2 [*<*(*» + 1 ) ~ m L 



\Y 



(19) 



The behavior of I m , a (B) between these two limits is 
shown as a function of magnetic field for 133 Cs in Figure 
|4| For positive m/,a the coupling increases monotonically 
before leveling off to the value (19), while for negative 



mj ]fl it increases with B, peaks, and then declines to the 
same value. At low fields, the coupling is approximately 
proportional to B, so that the resonance width is pro- 
portional to B 2 in this region. The range over which this 
behavior occurs is system-dependent; the coupling ele- 
ments for lighter alkali metals level off at smaller B than 
for Cs. 

The factor l/d/j, Tes in Eq. (16) produces wider reso- 



nances when the difference in slope between the bound 
and continuum states at J5 rcs is small. Particularly shal- 
low crossings and wide resonances can occur when there 
is a "double crossing" involving a bound state that just 
dips below the threshold (as a function of B) before ris- 
ing above it again. The magnetic fields at which this can 
occur are discussed in Section [Til CI below. 

The bound and continuum functions, ip n (R) and 
ipk(R), are eigenfunctions that correspond to different 
eigenvalues of the 1-dimensional radial hamiltonian Q9J) . 
They are thus orthogonal to one another, and the matrix 



FIG. 5. (color online). The functions contributing to the 
integral I n k of Eq. ( |18[ ) for CsYb. (a) Wavefunction for the 
n — —4 bound state, xp n =-i{R), which is the shallowest bound 
state for which crossings exist, (b) wavefunction for the low- 
energy continuum state, ipk(R)- (c) A£cs(-R). (d) the inte- 
grand of the matrix element, ij) n ij)uAC,. (e) The partial integral 

f H i> n (R')MR')M(R')dR'. 



element I n k of Eq. (18) is non-zero only because of the 
-R-dependence of A£jR) . 

Figure [5] shows how the integral /„& develops as a func- 
tion of R in a typical case. The upper three panels show 
ip n (R), ^pk(R) and A( a (R). Figure [5JJd) shows the inte- 
grand of Eq. ( 18 1, which is the product of the three. For 
weakly bound states, the bound and continuum func- 
tions remain almost in phase with one another across the 
width of the potential well, so that their product always 
maintains the same sign. The integral thus accumulates 
monotonically as shown in the bottom panel. Its value 
depends principally on A( a (R) between the inner turn- 
ing point and the potential minimum. Deeply bound 
states lose phase with the continuum at shorter ranges; 
in principle this produces some cancelation that reduces 



the value of the integral, but the effect of this is small for 
the near-dissociation levels considered here. 

Further insight may be gained by considering the 
integral I n k semiclassically. In the WKB (Wentzel- 
Kramers-Brillouin) approximation, the bound and con- 
tinuum wavcfunctions both oscillate with amplitudes 
proportional to k(R)" 1 ^ 2 in the classically allowed re- 
gion, where k(R) = [2fi(E - V(i?))/?i 2 ] 1/2 - For very 
weakly bound states and low collision energies, E may 
be neglected, so 



/>oo 

J„ fc cx / k(R)' l AC(R) dR, 



(20) 



where r m is the inner classical turning point at E = 0. 
This structure is clearly visible in Figure p[d) . 

Near threshold, the WKB approximation gives an in- 
correct ratio between the short-range and long-range am- 
plitudes of a scattering wavefunction. Quantum Defect 
Theory (QDT) [30] corrects for this using an energy- 
dependent function C(fc), which is 1 far from threshold 
but is given by 



C(k)- 2 = ka 



^-°f) 



(21) 



at limitingly low energy |3|. The correction amplifies 
the short-range wavefunction by a factor C(k) , which 
has a minimum value of (ka) 1 ' 2 when ab g = a but is 
approximately (fc/a) 1 ' 2 ab g when |a Dg | 3> a. 

Combining all these effects gives a semiclassical expres- 
sion for the Golden Rule width, 



(i a 

fi Obs 



gbgV 

a ) 



[I mf jB)}< 
NS^ rcs 



k(R)- 1 AC(R)dR 



(22) 



where N is the normalization integral for the WKB 
bound-state wavefunction, 



N = 



1 



fc(/?) _1 rfi?, 



(23) 



Rin 



which is taken between the classical turning points R ul 
and i?out at energy E n . Eq. (22) completely avoids the 



calculation of any quantal wavefunctions and gives results 
within 2% of the quantal Golden Rule width (16). 



The semiclassical approach may be taken one step fur- 
ther, with a small approximation. For a near-dissociation 
vibrational state with an interaction potential that varies 
as —CjR~i at long range, Le Roy and Bernstein [3T] have 
shown that the integral ( 23 ) is 



Rin 



k{R)~ l dR; 



irh 2 
~2/7 



C 



Ui 



r i 



\En\ 



j + 2 
2j 



(24) 



where T(x) is the gamma function. For the present case, 
with i?~ 6 , A is thus proportional to |_E„| 2 / 3 . Deeper 
bound states thus produce broader resonances, though 
generally at higher magnetic field. For the bound states 
of interest here, Eq. (24) is accurate to within 6%. 



As described below, different isotopes of Yb offer dif- 
ferent values of the scattering length <2b g . Eq. (22) shows 
that large values of A may occur when |obg] is either 
very large or very small: A is directly proportional to 
ab g when |obg| 3> ci, and inversely proportional to ab g 
when |db g | *C a. 



Overall the Golden Rule approximation ( 22 ) produces 



resonance widths that agree within 2% with those from 
full coupled-channel calculations. It also produces impor- 
tant insights into the origins of the widths, and makes it 
much easier to select systems and isotopic combinations 
with experimentally desirable properties. 



E. Sensitivity to the interaction potential 

The Feshbach resonance positions and widths are 
strongly dependent on the s-wave scattering length of 
the system. The background scattering length ab g , the 
binding energies of high- lying vibrational levels E n , and 
the non-integer quantum number at dissociation «d can 
all be related to a semiclassical phase integral <&(E), 



$(£) 



k(R)dR. 



(25) 



For a potential with long-range behavior V(R) 
—CqR^ 6 , the scattering length is 



a hg = a 1 - tan ($(0) 



(26) 



where a is the mean scattering length of Gribakin and 
Flambaum [52 , which is proportional to ([iCq) 1 / 4 . Val- 
ues for a for all the alkali metals with Yb atoms are given 
for representative isotopes in Table |HI| The non-integer 
quantum number at dissociation is 



,GF 



$(0) 



(27) 



where the superscript GF distinguishes the Gribakin- 
Flambaum value from the (less accurate) first-order 
WKB value (see section IIIC). It should be noted that 



ab g is a single- valued function of the fractional part of 



& 



and is independent of its integer part. 



Potential energy curves from electronic structure cal- 
culations for heavy molecules are typically accurate to 
at best a few percent. For curves that support 35 to 
70 bound states, such as those for the systems consid- 
ered here, this uncertainty is enough to span more than 
1 in vr>. It is thus not possible to predict ab g for these 
systems from electronic structure calculations alone. An 
experimental measurement is essential to limit the possi- 
ble range of Ob g - 





168 Yb 


i76 Yb 


6 Li 


36.29 


36.31 


7 Li 


37.66 


37.68 


23 Na 


50.50 


50.57 


40 K 


63.10 


63.24 


87 Rb 


74.52 


74.82 


133 Cs 


83.05 


83.48 



TABLE III. Mean scattering lengths a (in bohr) for the Alk- 
Yb systems. 



u D ( 172 Yb) A?; D (Yb) 


6 Li 


23 0.02 


7 Li 


25 0.02 


23 Na 


35 0.10 


40 K 


45 0.20 


87 Rb 


62 0.49 


133 Cs 


69 0.70 



TABLE IV. The integer part of vd for Alk-Yb systems, based 
on the potential curves from CCSD(T) calculations, together 
with the amount Ai>d by which vo may be tuned by varying 
the isotope of Yb. Note that the number of bound states is 

Vu + I- 



If the uncertainty in vu is much greater than 1 and 
we assume that the possible values of $(0) (and hence 
vd) are uniformly distributed over such a range of uncer- 
tainty in V(r), we find from Eq. (26) that there is a 50% 
probability that a^, g is in the range [0, 2a], and a 70.5% 
probability that it is in the range [—a, 3a]. 

Different isotopologues of the same molecule have dif- 
ferent reduced masses (i. Since k(R) is proportional to 



/i 1 / 2 , changing between different isotopes of Yb alters 
$(0), and hence vn and <2b g , in a very well-defined way, 
which depends only weakly on the potential well depth. 
For the case of LiYb, changing the heavy-atom isotope 
has very little effect on the reduced mass and therefore 
on ab g . For the heavier alkalis, by contrast, changing the 
Yb isotope allows the scattering length to be tuned over 
a wide range. Table |TV| summarizes the number of bound 
states and the amount by which it may be tuned for all 
the alkali-metal + Yb systems. 



III. RESULTS AND DISCUSSION 

We have previously calculated resonance positions and 
widths for the LiYb systems [21] , using estimates of a^g 
obtained from thermalization measurements for 6 Li 174 Yb 
[201 [22] . In the following subsections, we present calcula- 
tions of resonance positions and widths for two cases rep- 
resentative of the heavier alkali metals: RbYb, where the 
scattering lengths are approximately known, and CsYb, 



where the scattering lengths have yet to be measured. 



A. RbYb 

Interactions of RbYb mixtures have been studied by 
Gorlitz and coworkers [53"rf56"] . Baumer et al. [5J3 [M] 
measured thermalization rates and density profiles for 
mixtures of Rb with a variety of Yb isotopes, and in- 
terpreted the results in terms of background scattering 
lengths. In particular, 87 Rb 174 Yb was found to have an 
extremely large scattering length, which produced phase 
separation of the atomic clouds, while 87 Rb 170 Yb was 
found to have an extremely small one. Munchow et 
al. [55[ [55] measured 2-photon photoassociation spectra 
of high-lying vibrational states of the electronic ground 
state: for 87 Rb 176 Yb, 6 states were observed with bind- 
ing energies between about 300 MHz and 60 GHz [551156] . 
whereas for each of 170 Yb, 172 Yb and 174 Yb, two states 
were observed with binding energies between 100 and 
1500 MHz. Munchow [56] fitted the binding energies 
to a Lennard- Jones potential model and inferred from 
the mass scaling that the potential supports about 66 
bound states for 87 Rb 174 Yb and 87 Rb 176 Yb, with one 
fewer state for lighter Yb isotopes. The presence of a 
bound state very close to dissociation in 87 Rb 174 Yb pro- 
duces its large positive scattering length. 

The Lennard-Jones potential reproduces the experi- 
mental spectra satisfactorily, but the mass scaling de- 
termines only the number of bound states and there is 
no reason to expect the potential to have the correct 
well depth, equilibrium distance, or inner turning point. 
These features are however important in the calculation 
of resonance widths. We have therefore refitted the bind- 
ing energies measured by Munchow [56j . together with 
the scattering length for 87 Rb 170 Yb, to obtain a new po- 
tential curve based on our CCSD(T) results described 
above. Our best fit was obtained by multiplying the 
CCSD(T) potential by a scaling factor A sc] = 1.09581 
and adjusting Cq to 2874.7 E^a®, producing a potential 
that supports 66 bound states for 87 Rb 176 Yb. We also 
introduced long-range C% and Cio coefficients related to 
Cq by a ratio 7 = Cg/Cg = Cio/Cg, with an optimum 
value 7 = 267.5 bohr 2 for the potential above. Since 
the resulting long-range potential is valid to shorter dis- 
tances than the pure CqR~ & potential used for the other 
systems, the switching function [43] was applied between 
20 and 30 bohr in this case. It should be noted that 
adequate fits could also be obtained with one or two ad- 
ditional (or fewer) bound states: increasing A sc i by 0.036 
and Cg by 37 E^a® produces a potential with one extra 
bound state at the bottom of the well but the high-lying 
states almost unchanged. 

We have carried out coupled-channel calculations for 
the RbYb systems using the fitted potential with 66 
bound states. For the fermionic isotopes 171 Yb and 
173 Yb, we neglected couplings due to A£ Y b- The cross- 
ings responsible for the resonances for 87 RbYb are shown 
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FIG. 6. (color online). Resonance crossings for 87 Rb with the 
stable isotopes of Yb, demonstrating the mass-scaling effect. 
The threshold levels for rrif^h = — 1,0, +1 sublevels of the 
/ = 1 manifold are shown as dotted black lines, while the 
molecular bound state m^rtb = — 1,0,-1-1 sublevels for the 
/ = 2 manifolds are shown in different colors for the different 
isotopes of Yb. The bound states for m/,Rb = ~2 and +2 
are not shown. The highest bound state shown here is the 
n = —4 vibrational state for the combinations of 87 Rb with 
168 Yb^ 173 Yb, and n = -5 for 87 Rb 174 Yb and 87 Rb 176 Yb. 
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FIG. 7. (color online). Resonance crossings for Rb with 
the stable isotopes of Yb. The threshold levels for m/,Rb = 
—2, —1, 0, +1, +2 sublevels of the / = 2 manifold are shown as 
dotted black lines, while the corresponding molecular bound- 
state sublevels for the / = 3 manifolds are shown in different 
colors for the different isotopes of Yb. The bound states for 
Tnf,Kb = —3 and +3 are not shown. The highest bound state 
shown here is the n = —3 vibrational state. 



87-168 


-1 


3314 


4.6 


39 


1.4 


8.4 


x 10" 4 







1705 


2.0 


39 


1.1 


3.0 


x 10~ 4 




1 


877 


0.3 


39 


0.4 


5.9 


x 10" 5 




1 


4466 


2.4 


39 


0.5 


5.9 


x 10" 4 


87-170 


-1 


3723 


-31.1 


-11 


8.3 


1.8 


x 10~ 3 







2189 


-16.6 


-11 


7.5 


8.2 


x 10" 4 




1 


1287 


-3.8 


-11 


2.9 


2.1 


x 10" 4 




1 


4965 


-17.9 


-11 


3.6 


1.3 


x 10 -3 


87-171 


-1 


3923 


-10.9 


-58 


2.7 


3.4 


x 10 -3 







2415 


-6.2 


-58 


2.5 


1.7 


x 10" 4 




1 


1487 


-1.6 


-58 


1.0 


4.8 


x 10" 4 


87-172 


-1 


4122 


-10.4 


-156 


2.5 


8.9 


x 10~ 3 







2636 


-6.2 


-156 


2.3 


4.8 


x 10" 3 




1 


1686 


-1.7 


-155 


1.0 


1.5 


x 10" 3 



87-173 


-1 4320 


-20.7 


-576 


4.7 


6.7 x 10" 2 




2852 


-13.0 


-571 


4.5 


3.8 x 10" 2 




1 1883 


-3.9 


-569 


2.0 


1.3 x 10~ 2 


87-174 


-1 4517 


23.9 


991 


5.2 


1.4 x 10" 1 




3066 


16.0 


1000 


5.2 


8.5 x 10" 2 




1 2081 


5.2 


1005 


2.5 


3.0 x 10~~ 2 



87-176 



1 4912 


3.5 


224 


0.7 


4.7 x 10" 3 


3488 


2.5 


224 


0.7 


3.2 x 10" 3 


1 2476 


0.9 


224 


0.4 


1.2 x 10~ 3 



TABLE V. Predicted positions and widths for resonances with 
Am/, Rb = for 87 RbYb systems at fields B res < 5000 G. 



in Figure [6] and the resonance positions and widths are 
given in Table |V| for all resonances located below 5000 G. 
The corresponding results for 85 RbYb are given in Figure 
[7] and Table l*VTl for resonances located below 1500 G. A 
full listing of all resonances below 10000 G is provided 
as Supplemental Material [57]. The resonance positions 
are generally within about 50 G of those obtained by 
Miinchow [SS] with a Lennard- Jones model of the poten- 
tial. 

The pattern of widths for 87 RbYb closely follows ex- 
pectations from Eq. (I22J). Only 87 Rb 168 Yb has a res- 



onance below 1000 G, and that has a very low width 
(300 /J.G), in part because of dropoff in I m Rb (B) at low 
fields. Nevertheless, resonances with calculated widths 
as narrow as 0.2 /iG have been observed as 3-body loss 
features in Na [58 , and resonances a few mG wide have 
been observed in LiNa at fields as high as 2050 G (59) . 
87 Rb 170 Yb has particularly large widths as measured by 
A (up to 30 mG), but this is simply because Oh g is small 
in this case: the quantity ah g A, which is a better measure 
of the suitability of a resonance for magnetoassociation 
[60l l6T] , is not particularly large for this isotopologue. By 
contrast, 87 Rb 173 Yb and 87 Rb 174 Yb, which both have 
| ether | 3> a, have resonances up to 25 mG wide. Ex- 
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Rb-Yb 


Wl/,Rb 


-Ores 


A 


«bg 


|A| 

Srcs 




-Sres 






(G) 


(mG) (bohr) (ppm) 






85-168 


1 


1350 


0.17 


219 


0.12 


2.1 


x 10 -4 




2 


1100 


0.066 


219 


0.06 


8.8 


x 10~ 5 


85-170 


2 


1348 


0.048 


137 


0.04 


4.1 


x 10" 5 


85-171 


-2 


526 


-0.10 


116 


0.18 


1.7 


x 10~ 5 




-2 


918 


0.30 


116 


0.32 


5.4 


x 10~ 4 




2 


1475 


0.048 


116 


0.03 


4.1 


x 10~ 5 


85-172 


-2 


340 


-0.019 


99 


0.06 


5.4 


x 10~ 6 




-2 


1104 


0.21 


99 


0.19 


6.1 


x 10~ 5 


85-173 


-2 


206 


-0.0055 


84 


0.03 


1.5 


x 10~ 6 




-2 


1238 


0.21 


84 


0.17 


6.4 


x 10 -5 


85-174 


-2 


90 


-0.0009 


70 


0.01 


2.9 


x 10~ 7 




-2 


1354 


0.24 


69 


0.18 


7.0 


x 10~ 5 




-1 


270 


-0.11 


70 


0.39 


4.4 


x 10~ 6 




-1 


452 


0.30 


69 


0.18 


1.2 


x 10~ 5 


85-176 


-1 


925 


0.43 


39 


0.46 


5.5 


x 10~ 5 







434 


0.14 


39 


0.32 


1.4 


x 10~ 5 




1 


203 


0.021 


39 


0.10 


2.7 


x 10~ 6 




2 


120 


0.0033 


39 


0.03 


5.7 


x 10 -7 



TABLE VI. Predicted positions and widths for resonances 
with Am/.Rb = for 85 RbYb systems at fields B TCB < 1500 G. 



perimentally, 87 Rb 174 Yb displays phase separation that 
will inhibit molecule formation even for low-temperature 
thermal clouds [53], but 87 Rb 173 Yb does not [54], and is 
a good candidate for magnetoassociation if the high fields 
in Table [V] can be achieved. 

For 85 RbYb, there arc no resonances with A/i? ros > 
10~ 7 . This arises mostly because of the lower hyper- 
fine coupling constant £ for 85 Rb, which both reduces 
the magnitude of A£(i?) and further reduces the widths 
through the factor of |£" n | 2 / 3 described following Eq. (24) 



above. However, there are several resonances predicted 
below 1500 G, as shown in Table |VI| and some of the 
broader ones (still below 1 mG width) may be suitable 
for molecule formation. In particular, our best-fit po- 
tential predicts a pair of resonances for TO/jRb = —1 
for 85 Rb 174 Yb, where the atomic and molecular states 
just intersect and undergo a double crossing as shown in 
Fig. [7] The precise positions and widths of these res- 
onances are very sensitive to the potential details, and 
indeed Miinchow's Lennard- Jones model predicted that 
the atomic and molecular states just miss each other in- 
stead of just crossing [56 . 

Tables [V] and VI include only resonances driven by 
A((R) for Rb, which conserve m/m,. If the Yb 
isotope has nuclear spin, as for fermionic 171 Yb and 
173 Yb, additional resonances can occur at crossings with 
Am /3b = ±1, driven by AC(-R) for Yb [24]. In particu- 
lar, 87 Rb 171 Yb has a lower-field and therefore potentially 
more accessible group of resonances near 1210 G, where 
the molecular states with rrif, Rb = +2 and m^Yb (not 



shown in Fig. K3J) cross the thresholds with ro/^b = +1 
and Wi.Yb — 1- 

All the resonances in Tables [V] and VI are strongly 
closed-channcl-dominatcd. This may be quantified 
using the dimensionless resonance parameter s rcs = 
(ab g /a)(SfiA/E), where E — h 2 /(2^ia 2 ). It may be seen 
that s res is never greater than 0.2, and approaches such 
values only when |ab g | is very large. In some cases s res 
can be less than 10 -6 . 

Molecule formation by magnetoassociation is usually 
carried out by preparing the atomic mixture close to a 
resonance, on the side where the atomic state lies below 
the molecular state, and then ramping the field over the 
resonance. However, for narrow resonances in Cs2 (a few 
mG wide, at low fields), Mark et al. [62] found it effective 
simply to hold the field on resonance for a few millisec- 
onds. Nevertheless, the most efficient molecule produc- 
tion occurs with a field ramp that is slow enough to cross 
the resonance adiabatically (60] |6TJ [63] . Small field in- 
homogeneities are not a big problem, as they will simply 
cause different parts of the cloud to cross the resonance 
at slightly different times. However, field noise is poten- 
tially a problem, particularly high-frequency noise that 
causes nonadiabatic crossings through the resonance. It 
will therefore be important to design a molecule creation 
experiment with very careful field control. In this con- 
text it is worth noting that Ziirn et al. [64j have recently 
carried out radiofrequency spectroscopy on Li 2 molecules 
at fields around 800 G with a field precision of ±1 mG, 
which is close to 1 part in 10 6 , while Heo et al. [TS] 
achieved molecule formation in 6 LiNa, using a resonance 
10 mG wide at 745 G, with active feedback stabilization 
of the current to achieve field noise less than 10 mG 1181. 



B. CsYb 

Cesium possesses several properties that make it fa- 
vorable compared to the other alkali-metal elements for 
magnetoassociation with Yb. It has the highest mass 
of the alkali metals, which leads to greater mass scaling 
through changing the isotope of the closed-shell atom. 
Its larger mass also provides a higher density of bound 
states near threshold and thus offers better chances of 
resonances at low magnetic field. Additionally, its rela- 
tively large nuclear spin allows larger off-diagonal i a ■ s 
elements. Finally, the effects of A£(i?) are larger for Cs 
than for most of the other alkali metals. 

The CsYb potential shown in Fig.[2]supports 70 bound 
states for all Yb isotopes, and has a background scatter- 
ing length ab g = —38 bohr for 133 Cs 174 Yb. However, the 
electronic structure calculations have a degree of inaccu- 
racy, and a plausible change of ±10% in the well depth 
would produce a change of ±3 in vu- Since the scattering 
length depends on the fractional part of vu , it cannot be 
predicted from these calculations. However, altering the 
Yb isotopic mass across its possible range from 168 to 176 
changes «d by about 0.70, so that a wide range of back- 
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FIG. 8. (color online). Resonance positions and widths for 
133 Cs 174 Yb, obtained from coupled-channel calculations, as a 
function of the background scattering length ab g . Top panel: 
Resonance position B ICS . Center panel: Width |A|, with an 
inset showing an expanded view for values of ab g near the 
mean scattering length, a. Bottom panel: A/B rcs in parts 
per million. 



as Ob g increases from —2000 to +2000 bohr. When a 
state becomes too shallow to cross the lower threshold 
at all, the corresponding resonance line either disappears 
through B = or (in the case of a double crossing, as 
for TO/,Cs — — 3) reaches a maximum Ob g where the two 
crossings coalesce. 

The middle panel of Figure M shows |A| as a function 
of abg- The spikes in |A| near ab g = occur because of 
the Obg in the denominator of Eq . (|22[ ). However, the 
strength of the peak in a(B) [Eq. (|12[)j is actually ab g A 
rather than A itself, so this situation does not offer par- 
ticular advantages for molecule formation. For |ab g ^$> a, 
the widths vary linearly with abg as described in Sec- 
tion |IID| A particularly interesting feature of this plot 
is the spike in the to/,Cs = — 3 widths near ab g = —100 
bohr, which is physically significant. As noted above, 
the bound states for this magnetic sublevel experience a 
double crossing with the lower threshold in this region; as 
the two crossings approach one another, 5/j, Tes decreases 
and A increases as given by Eq. ( 16 ). A similar spike oc- 
curs in the to/.cs = — 2 resonance widths near 167 bohr. 
The inset shows an expanded view of |A| for the range 
of abg 



from —a to 3a; as described in section HE there 



is about a 70% probability that ab g lies in this range for 
any particular isotope. 



C. Choosing promising systems 

The Fermi Golden Rule treatment developed above 
shows that the most important properties leading to large 
resonance widths are a large magnitude of the back- 
ground scattering length and the occurrence of "double 
crossings" where the bound and continuum states have 
similar (small) magnetic moments. It is instructive to 
consider the conditions where these two enhancements 
can occur together. 

At zero field, the hyperfme splitting of an alkali-metal 
atom in a 2 S state is E^ (0) = £ (i + i) . As a function of 
magnetic field, the splitting between two states with the 
same value of m,f (neglecting the nuclear Zeeman term) 
is 



ground scattering lengths will be accessible by varying 
the Yb isotope. We have therefore carried out calcula- 
tions for CsYb as a function of abg- 

The top panel of Figure M shows the resonance posi- 
tions and widths for 133 Cs 17 ^b as a function of ab g - The 
plot would be almost identical for any other Yb isotope 
(though different isotopes will have different scattering 
lengths). There are multiple resonances for each value 
of rntfisi which occur when bound states \a2,rrif ,csi n ) 
cross the scattering threshold \a±, m/ : cs)- As Ob g in- 
creases, the binding energies decrease and the most of 
the crossings (those with positive Sfi) shift to lower mag- 
netic fields. In the nif^cs — +3 case, the position of one 
resonance changes from B w 6000 G to B w 2000 G 



E h( (B) = 



E M (0? 



4m f E hi (0) 
2i + l 



g e fi B B + (g e fi B BY 



(28) 



For negative values of mj, this has a minimum value 



E 



close 



E hl (0) 1 



Am'i 
(2z + l) 2 



at a field 



B 



■Aosc _ -2m/E hf (0) _ -to/C 



(29) 



(30) 



(2i + l)a e ^ B 9el^B 
The first-order WKB quantisation formula, expressed in 
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terms of the phase integral of Eq. (251, is 



(31) 



Le Roy and Bernstein [ST] showed that this implies that, 
for a long-range potential V(R) = —CjR~i with j > 2, 
near-dissociation levels exist at energies 



E„ 



■[Bj(vE 



WKB 



V)] 



2j/(i-2) 



where 



H, 



~2/7 



(i-2) 



r i + 



c 



i/j 



(32) 



(33) 



However, Eqs. (31) and (32) do not take account of the 
Gribakin-Flambaum correction [55], which replaces the 
(v + i) in Eq. (31) with (v + | + e(w)), where e(v) is 



zero for deeply bound levels but is | at dissociation for 
a long-range R~ e potential. This correction may have a 
significant effect on the energy of the least-bound level 
[651 166) , which is responsible for the Feshbach resonances 
Li-Yb [21] but is small for the slightly deeper levels that 
are responsible for the resonances in the heavier Alk-Yb 
systems. As a result, the near-dissociation levels (except 
n = —1) actually occur at energies close to 



E v = —H, 



,/„.WKB_^3 



= -Hi 



v§ F -v + 



where 



,,WKB 



$(0) 1 



,GF 



(34) 



(35) 



when expressed in terms of Up from Eq. ( 27 ) 



Very large values of |ob g | correspond to near- integer 
values of v*g F , so the condition for a very large value of 
|etbg| to coexist with a "double crossing" near B m 
that the dimensionless quantity 



lose 



IS 



v 

171 f 



ghf(0) 

m 



i - 



(2z + l) 2 



1/3 



(36) 



should be approximately an integer. The quantity X m , 
may be interpreted as the vibrational quantum number 
(relative to threshold) that will just give a double cross- 
ing between a molecular state associated with the upper 
hyperfine level and an atomic state at the lower hyperfine 
threshold with the same to/. It depends strongly on the 
alkali-metal isotope through the nuclear spin and hyper- 
fine splitting, but is only very weakly dependent on the 

1/2 

Yb isotope chosen. It is proportional to C 6 (through 
He), but is otherwise completely independent of the in- 
teraction potential. Values of X mj slightly smaller than 
than an integer allow a very large value of |ab g | to coex- 
ist with a double crossing further from BJjJ* c , or a large 



m f B£° f SB (G) 



x„ 



ar, x (bohr) a mm (bohr) 



"Li 


-1/2 


27 


0.20 


86 


84 


7 Li 


-1 


143 


0.40 


51 


47 


23 Na 


-1 


316 


1.06 


297 


175 


39 K 


-1 


82 


0.90 


-122 


-302 


40 R 


-7/2 


357 


1.17 


173 


87 


40 R 


-5/2 


255 


1.29 


111 


87 


40 K 


-3/2 


153 


1.35 


94 


87 


40 R 


-1/2 


51 


1.38 


88 


87 


41 K 


-1 


45 


0.71 


14 


-1 


85 Rb 


-2 


722 


2.36 


111 


48 


85 Rb 


-1 


361 


2.56 


61 


48 


87 Rb 


-1 


1219 


3.32 


123 


78 


133 Cs 


-3 


2460 


3.95 


-399 


71 


133 Cs 


-2 


1640 


4.33 


133 


71 


133 Cs 


-1 


820 


4.50 


84 


71 



TABLE VII. The quantity X mf of Eq. (36 1, which needs to 



be close to an integer for double crossings to exist for large 
values of |<2bg|, together with the range of scattering lengths 
for which double crossings can exist for a pure R~ 6 potential. 
The values are almost independent of the Yb isotope. 



negative value of a^ g to coexist with a double crossing 
near B* sc . In general, large negative values of Ob g may 
be more favorable for molecule formation than large pos- 
itive ones, because negative values will not cause phase 
separation in condensates. 

Values of X m . for all the Alk-Yb systems are given 



in Table |VII| They may be converted into values of the 
scattering length that just cause double crossings (for a 
pure Cq/R 6 potential) using 



= a {1 - tan [it (X mf + |)] } 



^ ax and 

1 J 



(37) 



will 



mm ^ ^max 



Scattering lengths between a m emu u — u 

give rise to double crossings (where values of a mm > a 
are to be interpreted as allowing the scattering length to 
be decreased from a™ ax , through a pole and back down 
from +oo to a min ). However, only values close to a™ ax 

result in double crossings close to B* so , which are the 
ones with particularly large widths. Table |VII| includes 
values of a™ ax and a min for all the Alk-Yb systems. It 
immediately explains why 85 Rb 174 Yb, with a background 
scattering length ab g = 70 bohr that is reasonably close 
to a™ ax = 61 bohr, can have a double crossing near 
B ciosc = 361 G for m/ = -1. The fact that this occurs 
with Ob g slightly larger than a™ ax (rather than slightly 



smaller) reflects the approximations inherent in Eq. ( 37 1 : 
it applies only to a pure Cq/R 6 potential and only ap- 
proximately incorporates the Gribakin-Flambaum cor- 
rection. Table |VII| also explains why Figure [8] shows 
peaks in resonance widths for CsYb at moderately large 
negative <Zb g for to/ = —3 and for moderately large pos- 
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itive abg for m/ = — 2. 

In general terms Yb is a favorable atom because it of- 
fers a large number of isotopes that facilitate tuning the 
reduced mass and hence ab g . The heavier alkali met- 
als are more favorable than the light ones because their 
larger masses offer greater tunability by varying the Yb 
mass. The heavier alkali metals are also more favorable 
because the levels that offer crossings at moderate mag- 
netic fields have larger binding energies (\E n \ between 
-E™ and -Em(O)). CsYb appears to be particularly fa- 
vorable because the near-integer value of X_ 3 makes it 
possible for shallow double crossings to coexist with large 
values of the scattering length. 

As discussed above, the short-range amplitude of the 
bound-state wavefunction is proportional to I-EVJ 1 ' 3 . In 
addition, A£ a (i?) is very roughly proportional to £ a : for 
the Alk-Yb systems, Co/Ca is about 0.3 for Li and Na and 
between 0.16 and 0.20 for K, Rb and Cs. The integral 
I n k of Eq. ( 18 ) thus scales very roughly as (, a for reso- 
nances that occur at fields below B closc . This effect itself 
accounts for a factor of nearly 20 between the resonance 
widths for 87 Rb and 85 Rb. 



IV. CONCLUSION 

We have investigated Feshbach resonances in mix- 
tures of alkali-metal atoms with Yb, in order to identify 
promising systems for magnetoassociation to form ultra- 
cold molecules with both electric and magnetic dipole 
moments. The resonances in these systems arise when 
molecular states associated with the upper hyperfine level 
of the alkali-metal atom cross atomic thresholds associ- 
ated with the lower hyperfine level. They are due to 
coupling by the distance-dependence A£(i?) of the alkali- 
metal hyperfine coupling constant [03]. The widths of 
the resonances range from a few microgauss to around 
100 mG. 

We have calculated the potential energy curves and 
A((R) for Yb interacting with Na, K, Rb and Cs. We 
have carried out coupled-channel calculations of the reso- 
nance positions and widths for all isotopologues of RbYb 
and CsYb, and have also developed a perturbative model 
of the resonance widths that gives good agreement with 
the coupled-channel results. Key conclusions of the 
model are (i) that resonance widths depend strongly on 
the atomic hyperfine coupling constant £, with a gen- 
eral scaling as C 8 ^ 3 ; (ii) that resonance widths are gener- 



ally proportional to the background scattering length ab g 
when it is larger than the mean scattering length a; (iii) 
that resonance widths are proportional to B 2 in the low- 
field region where the atomic Zeeman effect is linear; (iv) 
that unusually wide resonances may occur when a molec- 
ular bound state only just crosses an atomic threshold 
as a function of B; (v) that, for the heavier alkali met- 
als, varying the Yb isotope gives access to a wide range 
of background scattering lengths and thus to a range of 
different resonance positions and properties. Selecting 
the best isotope is likely to be crucial to the success of 
molecule production experiments. 

Accurate predictions of resonance positions and widths 
for a given system require knowledge of the background 
scattering length, or equivalently of the binding energy 
of the least-bound vibrational state. This cannot be 
obtained reliably from electronic structure calculations 
alone, and requires an experimental measurement on at 
least one isotopologue. Once this is available, the poten- 
tial energy curves from electronic structure calculations 
are accurate enough to allow mass-scaling to obtain pre- 
dictions for all isotopologues. For RbYb, for which bind- 
ing energies have been measured by 2-photon photoasso- 
ciation spectroscopy [SS], we have adjusted our potential 
curve to reproduce the experimental results and used the 
result to calculate resonance positions and widths. We 
find that some isotopologues of 85 RbYb have resonances 
at fields below 1000 G, but these are all very narrow 
(< 0.5 mG). Isotopologues of 87 RbYb have considerably 
wider resonances (some up to 30 mG wide) , but the most 
promising resonances occur at fields above 2500 G. 

For CsYb, no measurements of background scattering 
lengths or binding energies are yet available. We have 
therefore calculated the resonance positions and widths 
as a function of scattering length. CsYb is a particularly 
favorable combination because shallow double crossings 
may occur for isotopologues with large Ob g , producing 
particularly broad resonances. The mapping from scat- 
tering length to positions and widths is almost indepen- 
dent of isotopologue, although the actual values of ab g 
will be strongly isotope-dependent. 
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